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ABSTRACT 

High-resolution near-infrared HST NICMOS (F160W, F222M) images and polariza- 
tion (2 /xm) observations were made of four bipolar proto-planetary nebulae (PPNs): 
IRAS 17150-3224, IRAS 17441-2411, IRAS 17245-3951, and IRAS 16594-4656. The 
first three of these are viewed nearly edge-on, and for the first time the central stars 
in them are seen. Color maps reveal a reddened torus between the bipolar lobes in the 
edge-on cases, with bluer lobes. The polarization values are high, with maximum val- 
ues ranging from 40 to 80%. The polarization patterns are basically centrosymmetric, 
with some deviations in the low polarization equatorial regions. For IRAS 17150—3224, 
circumstellar arcs are seen at 1.6 /im, along with a newly-discovered loop in the equa- 
torial region. Bright caps are seen at the end of the lobes, indicating that they are not 
open-ended. A distinct point-symmetric pattern is seen in the strengths of the polar- 
ization vectors, especially in IRAS 17150—3224. HST NICMOS observations provide a 
valuable complement to the WFPC2 visible images in deriving the basic structure of 
bipolar PPNs. 

Subject headings: circumstellar matter — planetary nebulae: general — stars: AGB 
and post- AGB 



1. INTRODUCTION 

As low- and intermediate-mass stars enter the asymptotic giant branch (AGB) of evolution, 
they enshroud themselves in dusty cocoons created by heavy mass loss in the form of stellar winds. 
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Such winds gradually deplete the hydrogen envelopes of the stars, exposing their hot cores and 

causing the stars to increase in temperature. As the envelopes continue to thin as the result of 
nuclear burning and mass loss, the stellar temperature eventually increases to the extent that sig- 
nificant ultraviolet photons are emitted. The onset of photoionization leads to strong line emission 
from the circumstellar material, producing a planetary nebula (PN). 

In the past 25 years, significant progress has been made in the understanding of both the 
evolution of the central star and the dynamical evolution of the nebula. It is now believed that, 
in most cases, the PN central star is powered by hydrogen shell burning, and that the dynamics 
and the morphology of the nebula are driven by an interacting winds process (cf. Kwok 2000). 
While the one-dimensional structures (including shells, crowns, haloes, etc.) can be understood by 
the coupled stellar and nebular evolution models (Schonberner & Steffen 2003), the origin of the 
variety of two-dimensional structures is still an unresolved mystery. The observed morphologies of 
PN range from round to elliptical to bipolar to point symmetric, whereas their progenitors, the 
AGB stars, show spherical symmetry in their envelopes (e.g., Neri et al. 1998; Sahai k. Bieging 
1993; Lucas et al. 1992). Although it is certain that the interacting winds process plays a role in 
the shaping of the nebula (Balick 1987), the exact details on how and when such shaping occurs 
remain unknown. 

The study of proto-planetary nebulae (PPNs), the transitional objects between the AGB and 
PN phases, is crucial to an understanding of the processes that shape PNs. PN morphologies 
must be created or enhanced during this transitional phase. PPNs have recently been imaged in 
visible light with the Hubble Space Telescope (HST) and found to display many of the same shapes 
as the PNs (Ueta et al. 1998). Many show bipolar reflection nebulae, some with concentric arcs 
and "searchlight beams" (Sahai et al. 1998a; Kwok et al. 1998; Su et al. 1998; Hrivnak et al. 1999, 
2001). In these, there appears to be a higher density of material in the equatorial region forming an 
optically-thick dust cocoon, with low-density lobes in the polar directions. The bipolarity observed 
in PPNs suggests that the shaping of PNs is established early in the transitional phase. 

In this paper we present a high spatial resolution, near-infrared HST imaging and polarization 
study of four PPNs, each of which has recently been studied with high spatial resolution visible 
imaging. Three of these display a dark lane in the visible images that obscures the central star. In 
the near-infrared, one can probe the structure of the circumstellar shell closer to the central star 
than can be done in visible light. One of our observational goals was to detect the central stars 
and determine their precise locations. With the polarization observations we wanted to investigate 
the scattering geometry, locate the position of the central star from the polarization vectors, learn 
about the properties of the dust grains, and see if there was evidence to support the model of a 
circumstellar torus in the systems. The only PPN for which a detailed high-resolution near-infrared 
imaging and polarization study has been published is AFGL 2688 (Egg Nebula), and in this case 
the central star was not seen in the near-infrared (Sahai et al. 1998b; Weintraub et al. 2000). As 
part of our program of high-resolution near-infrared imaging and polarimetry of PPNs, we have 
also presented preliminary results for three other PPNs (Sahai et al. 2000). With an increased 
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sample, we hope to better understand the properties of the circumstellar material of PPNs and the 
cause of the development of asymmetry in PNs. 

2. OBSERVATIONS AND DATA REDUCTION 

High-resolution near-infrared observations of these four PPNs were obtained in 1998 with the 
Near-Infrared Camera and Multiobject Spectrometer (NICMOS) on board the HST under General 
Observer program ID No. 7840 (PI: Kwok). NICMOS camera 2 (NIC2) was used, which has a 
very high resolution ((X'076 x (X'075 pixel" ^) and a 19^'5 x 19^'3 field of view (Thompson et al. 
1998). Observations were made with a broad-band H filter (F160W), a medium-band K filter 
(F222M), and three 2 ^m (1.9-2.1 /im ) polarizers (POLOL, POL120L, POL240L). The data were 
read out using the multiple non-destructive mode (MULTIACCUM) to obtain a high dynamic range 
without saturating the detector. An observing log, including total exposure times in the individual 
bandpasses, is given in Table 1. Each target was imaged three successive times in each filter at 
different positions on the NIC2 array, using a predefined spiral dither pattern. The data from 
the bad pixels and columns and the area blocked by the coronographic mask were compensated 
by the three-step dither mode. The data were processed with the standard NICMOS CALNICA 
pipeline program, which includes bias subtraction, dark subtraction, flat-field correction, cosmic- 
ray removal, and flux calibration. An additional effect that needed to be corrected is the so-called 
"pedestal effect" due to residual charge present in certain images (Thompson et al. 1998). Through 
use of the "unpedestal" program developed by Roeland van der Marel of the Space Telescope 
Science Institute, most of this effect was removed. 

The dithered images were shifted and coadded manually using median combining. In order to 
combine the dithered images accurately, we first subgridded each data set by a subsample factor 
of two so that the peak position of field stars can be measured more precisely. The offsets between 
dithered images were determined by comparing the peak positions of a number of field stars and 
obtaining an average offset. The precision in this method of registering the dithered images can be 
seen from the typical rms values of the offsets in the peak positions of the field stars in an image, 
which are 0.040 and 0.034 pixels in the X and Y directions, respectively. 

The distortion of the NICMOS dewar due to cryogen expansion affected its optical path. This 
caused the effective focal length, and hence the pixel scale, to change somewhat as a function of 
time during the initial active lifetime of the instrument. The plate scales were measured every two 
weeks during the active lifetime; therefore, wc adopted a plate scale for each observation using the 
measurement that was closest to the date of the observation. The NICMOS detector arrays were 
slightly tilted relative to the focal plane; therefore, when projected onto the sky, the pixels are 
actually rectangular, not square. The X and Y pixel scale differences need to be corrected to carry 
out precision astrometry or to accurately register NICMOS images to images taken with other 
instruments such as the Wide Field Planetary Camera 2 (WFPC2) on board HST. These scale 
differences must also be accounted for when analyzing the orientation of the polarization vectors. 
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The images were thus rescaled spatially when necessary. 

To remove (or at least reduce) extended point-spread function (PSF) features of the NICMOS, 
we have deconvolved the images using an iterative Poisson maximum-likelihood algorithm (task 
"lucy" in IRAF). The PSFs used in the deconvolution were generated by the TinyTim program 
v.4.4 (Krist & Hook 1997). The deconvolved images were only used for display purposes. All 
scientific measurements were done based on the original images. Specific reduction and analysis of 
the polarimetric data are discussed separately in Sec. 4. 

3. F160W AND F222M IMAGES 

Photometric measurements were made on the F160W and F222M images of the four PPNs 
using the calibration values listed in Table 2. The magnitudes determined from the FITS image 
headers (PHOTFNU and PHOTFLAM) are in the HST STMAG system. Because the NICMOS fil- 
ter bandpasses do not match those of ground-based filters, it is not straightforward to derive stan- 
dard infrared magnitudes. Since there are no direct spectrophotometric observations of Vega or 
other AOV standards through the NICMOS filter system, we used the flux values for a 0.0 mag- 
nitude star (Vega) as derived by the NICMOS Instrument Definition Team (IDT; see Table 2). 
Standard aperture photometry was carried out on the four objects using elliptical apertures that 
approximately matched the size of the nebulae. Apertures were chosen that were similar in size 
to the nebular halos measured on the visible WFPC2 images. The several bright field stars inside 
the apertures were excluded in calculating the total integrated counts. Table 3 lists the results 
of the photometric measurements and includes the aperture sizes. These magnitudes on the Vega 
system, mvega(P160W) and myega(F222M), are all close (< 0.2 mag) to the ground-based and 
K measurements, respectively, for these objects. The brightness of each of the central stars was 
also measured. We used a small aperture (2 x FWHM) and applied an aperture correction based 
upon the measurement of a bright star in the field of IRAS 17150—3224. These approximate values 
are also listed in Table 3; we estimate the accuracy to be ±0.3 to ±0.4 mag, the main uncertainty 
being the contribution of the nebular flux. 

The published WFPC2 F606W {Xe/f = 0.58 fim) images of IRAS 17150-3224 (Kwok et al. 
1998), IRAS 17441-2411 (Su et al. 1998), and IRAS 17245-3951 (Hrivnak et al. 1999) show a 
prominent bipolarity in their nebulae, with a dark lane separating two opposite lobes in each. 
Similar bipolarity is also seen in these new near-infrared images. This indicates that there exists an 
asymmetry in their circumstcllar dust shells, with stellar light escaping from lower density regions 
to give rise to this bipolar nebulosity. To investigate further the presence of an asymmetry around 
each of these PPNs, flux-ratio or color maps were computed. For these we included the previous 
WFPC2 F606W images. These were shifted and rotated in order to align with the F160W and 
F222M images. This process was done with the IRAF/IMMATCH package, which uses three field 
stars that appear in all of the images as the reference points to perform geometrical transformations. 
Since the resolutions at different wavelengths and with different filter widths are not the same. 
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the higher-resolution WFPC2 F606W and NIC2 F160W images were convolved with an extended 

Gaussian profile to match the resolution in the NIC2 F222M image before computing the flux-ratio 
maps. These images and flux-ratio maps are presented and discussed for each object in the following 
subsections. 



3.1. IRAS 17150-3224 

Figures la and lb display the deconvolved F160W and F222M images of IRAS 17150-3224 
(Cotton Candy Nebula). The object consists of two bipolar reflection lobes that have similar sizes 
in all three bandpasses (F606W, F160W, and F222M). Bright "caps" appear at the ends of the 
lobes. The presence of these sharp caps strongly suggests that there exists a definite boundary to 
the material distribution in the polar directions rather than open-ended cavities. The bipolar lobes 
look like limb-brightened bubbles in the near-infrared images. Thus the density inside the lobes 
must be very low; they are essentially closed-end cavities. 

Superimposed on the bipolar lobes and the surrounding halo is a series of concentric, circular 
arcs. Six arcs are detected in the F160W image, and these correspond to the inner six arcs previously 
detected in the WFPC2 F606W image. This correspondence is clearly seen in the radial intensity 
profiles shown in Figure 2. An average separation of Of! 51 was obtained between the six arcs in the 
F160W image; this value agrees with the separation measured in the F606W image using a similar 
method. These arcs appear to be portions of incomplete circles. 

The central star can clearly be seen in both near-infrared images, in contrast to the visible 
image in which the lobes were separated by a dark lane with the star obscured. The position of the 
central star is similar in both, with an average value of R.A.(J2000) = 17''18™19';78, decl.(J2000) 
= — 32°27'21'.'2. This agrees with the position of the crossing point of the searchlight beams seen 
in the F606W image, R.A.(J2000) = 17'^18"*19?77 and decl.(J2000) = -32°27'21'.'l (Kwok et al. 
1998), and confirms that the cavity which produces the beams is centered on the central star. 

The inner boundries of the lobes in the near-infrared and visible images appear not to be 
parallel or axially-symmetric but rather open widely to the SW side (the right side in Figure la 
and lb). The inner bounderies appear to be brighter than the caps or the rest of the lobes in 
these near-infrared images. A faint loop extending from the central region toward the ENE and 
WSW directions is found in each of the near-infrared images; this loop is not seen in the visible 
image. The loop has an approximately elliptical shape centered at the central star, with a size of 
l'.'llx0'.'59 (semi-major axisx semi-minor axis) and a position angle (P.A.) of 63°. Note that the 
major axis of the loop is not perpendicular to the symmetry axis of the bipolar lobes (P.A. = 112°). 
This kind of loop structure has also been seen in some young PPNs (e.g., Hubble 5) in visible light. 

Figures Ic and Id display the color or flux-ratio maps for F160W/F606W (logarithmic scale) 
and F222M/F160W (linear scale). The bipolar reflection lobes and the faint arcs have bluer colors 
than the equatorial regions and the central star. The flux ratios in different parts of the nebula 
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are summarized in Table 4. The bluer colors in the bipolar lobes and the arcs are consistent with 
the idea that these regions are viewed in light scattered off the dust grains. The fact that the SE 
cap is redder than the NW cap suggests that the SE lobe is pointed away from us; thus its light 
is further reddened by passing through more of the halo. This is consistent with the visible image 
in which the SE lobe is fainter than the NW lobe. The reflected light from the bipolar lobes is 
mostly singly scattered, as deduced from the high percentage polarization seen at 2 //m (see Section 
4.3). This condition can be used to set constraints on the radial extinction of starlight and the 
line-of-sight scattering optical depth as a function of radius from the central star. The scattering 
optical depth at any radius can be estimated by comparing the observed scattered light intensity 
to the incident stellar radiation (e.g., Sahai et al. 1999b). The observed brightness for the brighter 
NW lobe at 2.22 /xm is 3.5xl0~^^ ergs cm~^ arcsec"^ at the projected radius of 2'.'7. The 
incident stellar radiation at 2.22 /xm (2.85 xlO~^ ergs cm~^ /xm~^) can be estimated using the 
bolometric flux of 5.71 xlO~^ ergs cm~^ s~^ (Kwok et al. 1996) and T(,ff = 5500 K. Assuming a 
1/r^ density distribution, eqn (1) from Sahai et al. (1999b) gives a line-of-sight optical depth of 
Tios =2.2x10^^, which is about 200 times smaller than we observed (T2.2 = 4.2) based on the flux 
ratio of the whole nebula (see the calculation below). This suggests that the integrated light from 
the nebula is mostly reddened due to the radial extinction from the star to the point where it gets 
scattered, and then slightly blue-ed by the scattering process. 

The central star is extremely red, with approximate color indices of {V—H) > 9 and {H—K) 
^ 2.0. The spectral type of the central star, based upon observations of scattered light from the 
lobes, is G2 I (Hu et al. 1993). The intrinsic color for an unreddened G2 I star is {V — H)=1.5 and 
(H — K)=0.l; thus the central star is highly reddened. The central star shows a F222M/F160W flux 
ratio of 1.2, which is signiflcantly larger than that expected (0.46) from its G2 I spectral type and 
estimated interstellar extinction (Ay ~ 1.0 mag; Su et al. 2001). This additional reddening can be 
attributed to circumstellar material, particularly in the dense equatorial region. By comparing the 
expected value of the F222M/F160W flux ratio (0.36) from an unreddened G2 I star {T^ff = 5500K) 
to our observed value (1.2), (assuming the standard interstellar extinction curve (Cardelli et al. 
1989) for the nebular dust), we find an extinction optical depth at 2.2 iim of T2.2 = 1-8 and at 
1.6 //m of ri.6 = 3.0 toward the central star. However, the total observed flux from the object 
gives a F222M/F160W flux ratio of 6.3, which gives r2.2 = 4.2 and ri.g = 7.1. The ti,q/t2.2 ratio 
from the central star (1.6) is smaller than that derived from the whole object, suggesting that the 
extinction curve between 1.6 and 2.2 ^m for dust grains along the equatorial direction is flatter 
than the interstellar one. This flatter wavelength dependency on the extinction curve is likely due 
to larger grain sizes in the equatorial dense region than the typical value (0.1 ^m ) for interstellar 
grains. Other PPNs, like the Egg Nebula (Sahai et al. 1998a) and IRAS 04296-^3429 (Sahai 1999), 
also show evidence for large grains in their dense equatorial regions. 

The flux-ratio maps reveal non-axially symmetric patterns, with the region between the lobes 
appearing as a red wedge opening to the SW (right side of Figure Ic and Id). This suggests that 
the dense equatorial region (torus) is larger in vertical extent on this side. The F222M/F160W 
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fiux-ratio map shows a red source at the location of the central star, with additional red regions 

nearby on each side. To better illustrate the structures in the central part, three intensity cuts ((/.'S 
wide) through the central red source were made near the x-axis of the F222M/F160W flux-ratio 
map (Fig. 3). The nearby red regions appear as peaks in these cuts, with two pairs of peaks (radii 
~(X'5 and and perhaps a third (~1'.'6) seen, in addition to the central peak. The pair of 

peaks at radius ^1" is associated with the extended loop seen in the near-infrared images. The 
red color presumably arises from the scattering of the extremely reddened light of the central star. 
The possible outermost pair of peaks located at r~l'.'6 is seen primarily along one cut. These peaks 
arise from regions where the S/N is low in both the F160W and F222M images. It is possible that 
they represent a second, outer loop, but at present we consider the physical reality of the outermost 
peaks to be uncertain. The inner pair of peaks may represent a cross-section of the torus. However, 
they have different angular separations from the central star (~0'.'34 for NE and ~0'.'52 for SW). 
If the inner pair represents the dusty torus, the non-equal separations imply that either the star is 
not at the center of the torus or the torus is not axially symmetric. The former of these can arise 
if the post-AGB central star is in orbit with a binary companion. (Such a companion would likely 
be a main-sequence or white dwarf star and therefore much fainter and not visible.) The average 
separation of the inner peaks translates to 6.4x10^^ cm (D/1 kpc). 



3.2. IRAS 17441-2411 

Figures 4a and 4b show the deconvolved NIC2 F160W and F222M images of IRAS 17441-2411 
(Silkworm Nebula). This object showed a pronounced S-shaped nebula in the visible F606W image 
(Su et al. 1998), and this S-shaped morphology is seen in both of the near- infrared images. This 
morphology suggests point-symmetry in the nebula. The central star, which was not seen in the 
visible, is clearly resolved in both of the near-infrared images, with a position of R.A.(J2000) = 
17''47™13!51, decl.(J2000) = -24°12'5l'.'6. This position agrees with the center position of the 
two pairs of concentric arcs seen in the F606W image, R.A.(J2000) = 17'*47'"13?49, decl.(J2000) 
= — 24°12'51'.'0. The nebula has a similar size in all three bands; this is to be expected when one 
is seeing scattering from the dense-walled lobes. 

Figures 4c and 4d show the flux-ratio maps. The reflection lobes are bluer than the equatorial 

region, as is the case in IRAS 17150^3224. The southern lobe is redder and fainter in visible light 
than the northern lobe; both of these characteristics suggests that the southern lobe is tipped away. 
The detailed flux ratios in different regions are listed in Table 5. The central star is very bright in 
the near-infrared images, resulting in a complex PSF; hence the flux-ratio map of F222M/F160W 
is complicated and full of PSF features (wings and spikes). In an inset in Figure 4c, the intensity 
contours around the central star in the F160W/F606W flux-ratio map are shown. These con- 
tours appear to show an axisymmetric torus perpendicular to the lobes. Thus, similar to IRAS 
17150—3224, the color distribution in IRAS 17441—2411 supports the model of a circumstellar 
torus. The central star appears extremely red (presumably as a result of circumstellar extinction). 
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with approximate values of {V—H) > 9 and {H—K) ~ 2.5. 

3.3. IRAS 17245-3951 

Near-infrared images of IRAS 17245—3951 (Walnut Nebula) are shown in Figure 5. The central 

star, which was not seen in the visible, is clearly resolved in both of the near-infrared images, with 
a position of R.A.(J2000)=17''28™04!59, decl.(J2000)=-39°53'44'.'9. The nebula is smah and has 
a size in the F160W image similar to that seen in the F606W image. However, the nebula is not 
seen in the F222M image; this may be due to the fact that the image does not go as deep as the 
other two and partly perhaps to the relatively large effect of the complicated PSF of the bright 
central star in the F222M image. A jet-like structTirc seen in the southern lobe in the F606W image 
(Hrivnak et al. 1999) does not appear in the F160W image. 

Figure 5c shows the F160W/F606W flux-ratio map. The bluest colors are associated with 
the lobes and the reddest colors with a torus perpendicular to the lobes. The northern lobe is 
bluer than the southern lobe, again consistent with the visible image in which the northern lobe is 
brighter than the southern lobe. The central star is very red, with ( V—H) > 3.5 and (H—K) ~ 1.6. 

3.4. IRAS 16594-4656 

The visible F606W image of IRAS 16594-4656 (Water Lily Nebula; Hrivnak et al. (1999)) 
indicates that the bipolar nebula is oriented at an intermediate inclination, not nearly edge-on as 
in the other three cases. The central star is seen in the F606W image, and it has been assigned a 
spectral type of B7 (Van de Steene, Wood, k van Hoof 2000). It is reddened, with (V-H) ~ 4.2 
and {H—K) 0.3. No dark lane is visible. Near-infrared images of IRAS 16594—4656 are shown in 
Figure 6. The nebulosity in the F160W image bears a general resemblance to the one seen in the 
visible image, but does not show as much faint detail. In the F606W image, included for comparison 
(Fig. 6c), the main lobes are along the ESE and WNW directions with some extended structures 
("petals") along the NE and SW directions. The nebulosity seen in F160W is also elongated along 
the ESE and WNW directions, and the petals appear very faintly. Very little nebulosity is detected 
in the F222M image; this is due in part to the fact that the PSF of the bright central star masks, to 
some extent, the light of the nebula. In the F160W/F606W flux-ratio map, the faint petals stand 
out as having a bluer color than the rest of the nebula. 
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4. 2 POLARIZATION IMAGES 

4.1. Reduction and Analysis Procedures 

The NICMOS polarization data of the four PPNs were analyzed using the algorithm developed 
by Hincs et al. (2000) to derive the Stokes images {I,Q, U) and the maps of percentage polarizations 
{p) and position angles {9). Before constructing the final Stokes images and percentage polarization 
and position angle maps, there were several additional steps in the data reduction. Any sky 
background was removed, since this would produce a small DC offset that would introduce some 
artificial polarization. The sky values in each image were determined by examining several blank sky 
areas and an average sky value was subtracted from the image. The X and Y pixel scale differences 
were adjusted to produce square pixels; otherwise, a systematic offset in polarization position angle 
would have been introduced (Mines et al. 2000). In order to transform the data to square pixels, 
the flux of each data set was redistributed with a bicubic spline interpolation function using a task 
called "imlintran" in IRAF. Since the percentage polarization and position angle maps are formed 
by dividing one image by another image, one must be careful to avoid a run-away error, which can 
occur when dividing by a very small number with a relatively large uncertainty. Therefore, before 
calculating the result maps, the three data sets were rebinned by 2x2 pixels to smooth out small 
fluctuations. Then the threshold value for a real signal was defined as 10 a sky'-, any signal below 10 
Csfej/ was considered as noise with zero percentage polarization. 

The actual uncertainties in the polarization will be larger than the Poisson errors. Additional 
sources in the form of systematic errors include the background subtraction and the image align- 
ment. In order to reduce possible errors in the final percentage polarization and position angle 
maps, a threshold of S/N = 10 was also set in the final Stokes images. As the following estimations 
show, with this conservative threshold, the estimated uncertainties are not important. We first 
added 1 cTsky noise to the sky-subtracted images (POLOL, POL120L, POL240L) and found the per- 
centage polarization to differ only within ±3%. We also computed a new set of polarization maps 
using un-shifted data and found that the percentage polarization differs by less than 3% within the 
the S/N=10 threshold region, although it can increase up to 24% outside that boundry. According 
to the on-orbit calibration of the NICMOS polarimetric capability, a total instrumental polarization 
of < 1% was found (Hines et al. 2000). Since these four objects are all in the general direction of 
the center of the galaxy (/: 340° to 4°, h: —3.3° to -|-3.0°), a small amount of polarization (<5%) 
may be of interstellar origin. Therefore, we conclude that, after correcting for the projection effect, 
rebinning, and setting the above thresholds, the uncertainty in the polarization is < 5% and ae 
< 2° (the later taken from Hines et al. (2000)). 
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4.2. Polarization Structures 

The polarization pattern for a simple reflection nebula with a central illuminating source, such 
as a PPN, is expected to be centrosymmetric about the illuminating star. Even for a slightly 
more complicated system composed of a mixture of reflected central starlight and line emission 
from the extended nebulosity, sTich as is found in a young PN, one still expects to observe this 
centrosymmetric pattern. This is because the line-emission radiation will be unpolarized and its 
effect will only be to dilute the polarized scattered light (Scarrott &: Scarrott 1995). The most 
likely near-infrared emission line in the PPN phase is from H2, and the NIC2 polarizer bandpass 
may admit as much as 10 % of any H2 emission present at 2.12 ^um. Among these four PPNs, only 
IRAS 17150—3224 is known to have H2 emission. This dilution effect will be discussed in 4.3. 

Figure 7 shows the 2 nm polarization maps for the four PPNs. In this figure, the polarization 
vectors are plotted as straight lines whose lengths represent the strength of the polarization; these 
arc superimposed on the contour images of the total flux of the light through the three polarizers 
(Stokes / image). The distribution of the polarization vectors shows a general centrosymmetric 
pattern around the central source in each of these PPNs. This is consistent with what one would 
expect from a reflection nebulosity with a central illuminating source. Some departure from a 
centrosymmetry pattern is found in the equatorial regions. This is commonly seen in the equatorial 
regions of bipolar reflection nebulae (Kastner &; Weintraub 1995; Scarrott & Scarrott 1995; Sahai et 
al. 1998b). In these equatorial regions, the polarization is low (< 10 %), as will be discussed below. 
This suggests that multiple scattering has reduced the strength of the polarization and provides 
further evidence for the presence of a dense circumstellar torus in these PPNs. 

In addition to the centrosymmetry, the polarization vectors also show a strong point symmetry 
in the distribution of polarization strengths around the central star. This can be seen in plots of the 
polarization vectors within different ranges of percentage polarization. IRAS 17150—3224 shows 
remarkable point symmetry, as seen in Figure 8. The lobes and the two caps are highly polarized, 
with p > 60%. The maximum percentage polarization is found at the edge of the SE lobe (the 
fainter lobe in the optical image), with p 76%. In general, the light from the lobes is polarized 
with p > 50% (see the combination of Fig. 8a and 8b). Lower polarization is found along the edges 
of the lobes. The areas with p < 10% are located in the equatorial region. For IRAS 17441—2411 
(Figure 9), the reflection lobes typically possess p ~ 45% polarization. The maximum polarization 
is found at the bend of the S-shape in the southern lobe (the fainter lobe in the optical image), 
with p ~ 60%. Similar to IRAS 17150—3224, lower polarization is found along the edges of the 
lobes and the areas associated with the lowest polarization (p < 10%) are located in the equatorial 
region. 

For IRAS 17245-3951 (Figure 10), the higher polarization (42% > p > 30%) is found near 
the ends of the nebular lobes. Again, the low polarization {p < 10%) region is associated with the 
equatorial direction. For IRAS 16594—4656 (Figure 11), the higher polarization region (59% > p > 
40%) is found to be associated with the edges of the nebula and the maximum polarization is found 
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in the eastern lobe. The vectors also show a general centrosymmetric pattern. Since the nebula is 
oriented at an intermediate inclination (Su et al. 2001), the torus is no longer edge-on and the low 
polarization (p < 10%) region is distributed around the central star. 

4.3. Polarization Discussion 

In Figure 12 are shown the percentage polarization maps for the IRAS 17150—3224 and IRAS 
17441—2411. One can clearly see the highly polarized regions at the ends and the edges of the lobes 
in IRAS 17150-3224. Note that H2 emission is observed in IRAS 17150-3224 (Su et al. 2000). 
The strongest H2 emission is seen from the polar caps, which are among the regions with the 
largest polarization strength (~ 60—76 %). However, the total flux contributed by the II2 emission 
is low compared to the total flux through three polarizers (~2.3 %). The undiluted percentage 
polarization should be higher by only ~1 %. 

Comparison of our HST imaging polarimetry with that obtained from the ground shows the 
advantages of the high ffS'T spatial resolution when observing small nebulae such as PPNs. Gledhill 
et al. (2000, 2001) carried out ground-based near-infrared (J-band) imaging polarimetry for three 
of these PPNs. They found the polarization vectors distributed in an elliptical pattern aligned 
orthogonal to the polar axis of each nebula, rather than in the general centrosymmetric pattern 
found in these HST NICMOS observations. This can be quite easily understood as the result 
of the lower spatial resolution (by a factor of five) of their observations, which results in spatial 
smearing of the vectors (see Lucas & Roche (1998), Fig. 12). Similarly, their values (at 1.25 /xm) 
for the maximum observed polarization are less than these HST values (at 2.0 //m). This is also 
likely due to spatial smearing of the polarization values and to the greater effect of dilution on the 
polarization values by the unpolarized light of the central stars in their lower-resolution images. 
These effects will be greatest for the small nebula of IRAS 17245—3951. (This difference in the 
maximum observed polarization values is not due to the difference in wavelengths observed, since 
these authors also observed some PPNs in both J and K (2.2 ^um) and found the maximum observed 
polarization at K to be comparable to or less than that at J). 

The polarization maps depend on the dust grain properties, in additional to the density struc- 
ture of the nebulae. We therefore use the Monte Carlo multiple-scattering models developed for 
bipolar nebulae (Whitney & Hartmann 1993; Whitney et al. 1997; Lucas & Roche 1998) around 
young stellar objects (YSOs), since these are morphologically rather similar to PPNs for which 
such models are not available. With these we analyze our polarimetric images and infer rough 
constraints on the dust properties of our PPNs. 

First, one can get some constraints on the grain albedo. Lucas Sz Roche (1998) simulated K- 
band images for low albedo {lok = 0.22) and high albedo {ojk = 0.8) dust grains (see their Fig. 14). 
The low albedo image shows pinched-in, bipolar-like intensity contours in the equatorial plane for a 
near edge-on system, while the high albedo image shows more elliptical-like contours. Since IRAS 
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17150—3224 and IRAS 17441—2411 appear bipolar, low albedo grains are likely to be dominant in 
these two. Furthermore, a lower albedo in the K band is consistent with single scattering of the 
photons, since a high albedo would lead to multiple scattering and a lower value for the maximum 
observed percentage polarization. The values of the observed percentage polarization are very high 
for these PPNs. Both the high values of the observed percentage polarization and the pinched-in 
morphologies infer that low albedo grains are dominant in these PPNs. A similar result was also 
found in another PPN, Roberts 22 (Sahai et al. 1999b). 

The value of the maximum linear percentage polarization (pmax) is a function of the grain size 
and thus can be used to constrain the size of the grains present. The value of Pmax is the highest 
percentage level of polarization that can be produced by single scattering, and the maximum 
observed value of the percentage polarization provides a lower limit to Pmax- In Table 6 are listed 
the observed maximum values of the polarization (pmax{oi>s)) for these four PPNs and the Egg 
Nebula (for comparison). Lucas Sz Roche (1998) calculated the percentage polarization produced 
by scattering though 90° using particles of different radii a (see their Fig. 13). They found that for 
spherical silicate grains (with refractive index m = 1.7 — O.OSi), Pmax is 100 % for grains with radii 
a < 0.3 fim but then drops off rapidly for larger grains. (For a metallic dust mixture (m = 1.7 — 
0.6i) this size constraint is removed since they maintain high polarization for a wide range of sizes.) 
Thus, if we assume silicate grains (silicate absorption at 9.7 //m is observed in IRAS 17150—3224; 
Kwok et al. (1995)), then these high observed maximum polarization values indicate that the dust 
must be dominated by small grains in the scattering envelope (the bipolar lobes and halo). Table 
6 also includes the location in the bipolar nebulae where the observed maximum polarization is 
found and the nebular size in visible light. Except for IRAS 17245—3951, the observed maximum 
polarization is found associated with the fainter lobe for each of these bipolar nebulae. In an 
axisymmetric distribution of dust illuminated by a central source and inclined to our line of sight, 
if the fainter lobe is more polarized than the brighter lobe, the most likely explanation is that 
the grains are forward scattering. This is because with forward scattering, the scattered flux is 
weighted toward small scattered angles (Whitney & Hartmann 1993). Thus the comparison of 
these near-infrared images and polarization maps with models developed for YSOs suggests that 
the dust in the scattering envelope is dominated by small, forward-scattering grains (radii a < 0.3 
/Ltm for silicates) with low albedos. 



5. SUMMARY AND CONCLUSIONS 

A high-resolution near-infrared imaging and polarization study of four bipolar PPNs has been 
carried out using NICMOS on the HST. In the three nearly edge-on nebulae, the central stars are 
seen for the first time in these near-infrared images. (This is in contrast to the Egg Nebula, which 
appears similar in visible light but in which the central star is not seen in the near-infrared.) All 
four show a general ccntrosymmetric pattern in their polarization vectors, indicating that one is 
seeing light from the central star that has been singly-scattered from the dust grains. The nebular 
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size is the same in the near-infrared images as in the previous visible HST images, consistent with 
scattering from dense-walled lobes. Color maps show that the light from the lobes is bluer than 
from the central region, and the polarization images show that the lobes are highly polarized (>40% 
for IRAS 17150-3224 and IRAS 17441-2411). These all indicate that the lobes are seen almost 
entirely in scattered light and that any emission sources in the lobes are faint. The dark region 
between the lobes is very red in the color maps and the polarization there is weak and deviates 
from the centrosymmetric pattern; these are consistent with the presence of a dense torus that lies 
perpendicular to the bipolar lobes and causes multiple scattering in the emergent light. 

The limb-brightening of the lobes seen in IRAS 17150—3224 suggests that they possess low 
density and are viewed primarily in light scattered from the cavity walls. Furthermore, the cap 
structures indicate that these cavities are not open-ended cones. An extended loop was discovered 
along the equatorial direction in this object. 

In addition to the centrosymmetric polarization patterns, point-symmetric patterns are also 
seen among the polarization vectors for these four PPNs. This is particularly striking for IRAS 
17150—3224, the largest of the edge-on nebulae, and was seen previously in the Egg Nebula (Wein- 
traub et al. 2000). Such point-symmetric structures resemble the BRET (bipolar, rotating, episodic 
jets) phenomenon observed in PNs (Lopez et al. 1995), and suggest that the mechanism responsible 
for producing such structures is already operating in PPNs. Although it is unclear at this time 
whether the BRET phenomenon is related to the multipolar lobes observed in PN 2440 (Lopez et 
al. 1998) and in low-excitation (and young) PNs (Sahai &; Trauger 1998), our present data suggest 
that colliminated bipolar outflows commence early in the post-AGB evolution. 

The results in this paper increase the number of detailed high-resolution near-infrared imaging 
and polarization studies carried out with the HST from one (the Egg Nebula, Sahai et al. 1998b) to 
five. Our observations confirm the basic model of PPNs in which an equatorial region of enhanced 
density (torus) obscures the central star in visible light and in which light escaping in the polar 

directions scatters from the walls of polar cavities carved out by a fast wind. These results agree 
with other recent observations that show that the morphological shaping of PN occurs early, before 
the onset of photoionization. This model can be further tested by comparing the high-resolution 
mid-infraxed images of PN and PPNs that are starting to be obtained with the new 8—10 m 
telescopes. 
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Table 1. HST NICMOS Observing Log 



IRAS ID 


Observation 
Date 




Total Exposure Times (s) 




F160W 


F222M 


POLOL 


POL120L 


POL240L 


16594-4656 


1998 May 2 


256 


288 


288 


288 


288 


17150-3224 


1998 Aug 16 


448 


416 


448 


448 


448 


17245-3951 


1998 Apr 30 


256 


272 


272 


272 


272 


17441-2411 


1998 Mar 10 


416 


416 


416 


416 


416 



Table 2. HST NICMOS Photometric Calibration Values ^ 



Filter 




AA 


PHOTFNU 


PHOTFLAM 


Fj/ for 0.0 mag 


PSF FWHM 




{llTQ. ) 


(/im ) 


(Jy/(ADU s-i)) 


(ergs cm~^ s"^ /xm 


(Jy) 


(") 










/(ADU s-i) 






F160W 


1.593 


0.277 


2.190E-06 


2.589E-15 


1083 


0.16 


F222M 


2.216 


0.120 


5.487E-06 


3.352E-15 


668 


0.20 



'^The photometric calibration was done by the NICMOS IDT (M. Rieke, private communication). 
^Based upon a Gaussian fit to the PSF. 
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Table 3. HST NICMOS Photometry Results 



IRAS ID 




Entire PPN 






Central Star 




Size^ 


FIGOW 


F222M 


F160W 


F222M 




(") 




111 \ 


m.v;/'' 


iin ■(,/</ 


111 \ 


lUW.y/ 


16594-4656 


6.15x4.28 


12.57 


8.94 


13.16 


8.28 


8.9 


8.5 


17150-3224 


8.97x4.83 


13.67 


10.04 


14.13 


9.30 


13.2 


11.5 


17245-3951 


3.07x2.43 


14.11 


10.48 


14.63 


9.76 


11.7 


10.7 


17441-2411 


5.85x4.12 


13.70 


10.07 


14.20 


9.33 


13.2 


11.1 



^Size of the elliptical aperture: major axis x minor axis. 
^Calibrated magnitude on the HST STMAG system. 
^Standardized magnitude with respect to Vega. 



Table 4. Flux Ratios of the Various Components of IRAS 17150-3224 



Flux Ratio 


Lobes 


Arcs 


Halo 


Caps 


Central Star 


Torus 










NW 


SE 




NE 


SW 


F160W/F606W 


0.3~0.6 


0.5 


2.0 


1.3 


3.2 


>65 


24 


33 


F222M/F160W 


0.3 






0.3 


0.5 


1.2 


1.2 


1.3 



Table 5. Flux Ratios of the Various Components of IRAS 17441-2411 



Flux ratio 


N Lobe 


S Lobe 


Halo 


Central Star 


F160W/F606W 


1.8 


2.3 


1.8 


>60 


F222M/F160W 


0.36 


0.42 


0.42 


2.0 
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Table 6. Maximum Observed Percentage Polarization Values 





IRAS 


IRAS 


IRAS 


IRAS 


Egg Nebula 




17150-3224 


17441-2411 


16594-4656 


17245-3951 




Pmax{ohs) 


76 % 


60 % 


60 % 


42 % 


96 % 


Location 


fainter lobe 


fainter lobe 


fainter lobe 


brighter lobe 


fainter lobe 


Size (") 


9.3x3.4 


5.4x2.7 


8.4x4.7 


3.0x2.0 


12x4 
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Fig. 1.— Near-infrared images of IRAS 17150-3224: (a) F160W, (b) F222M, (c) flux-ratio map 
of F160W/F606W, and (d) flux-ratio map of F222M/F160W. The first three are displayed with a 
logarithmic scale and the fourth with a linear scale. In the flux-ratio maps, the redder regions are 
darker. The intensity range in (a) is 10~^'^ to 10~^ ergs cm~^ s~^ /xm ~^ arcsec"^ and in (b) is 
10~^^ to 10^^° ergs cm^^ s^^ fim arcscc^^. An inset in the F160W image uses intensity contours 
to show the region around the central star. 

Fig. 2.— Radial intensity profiles of the arcs in both the WFPC2 F606W and NIC2 F160W images 
of IRAS 17150—3224. The first six arcs (A— F) seen in the visible WFPC2 image are also clearly 
seen in the near-infrared NIC2 image. The arc F (peak value) has a 3 sigma-detection in the 
de-convloved H image. 

Fig. 3.— Intensity profiles along three cuts in the F222M/F160W color map of IRAS 17150-3224. 

Left: enlarged image of the central region with the three cuts indicated. Right: intensity profiles 
along the three cuts. The central star is at zero radius: positive radii represent the right side of the 
torus (SW) and negative radii represent the left side of the torus (NE). The position angles (P. A.) 
are measured with respect to the horizontal direction. 

Fig. 4.— Near-infrared images of IRAS 17441-2411: (a) F160W, (b) F222M, (c) flux-ratio map of 
F160W/F606W, and (d) flux-ratio map of F222M/F160W. All are displayed on a logarithmic scale. 
In the flux-ratio maps, the redder regions are darker. The intensity range in (a) is 10^^^ to 3x10^^ 
ergs cm^^ s~^ /im ~^ arcsec"^ and in (b) is 3x10^^^ to 3x10"^*^ ergs cm^^-s^^ //m arcsec^^. An 
enlargement of the equatorial region of the F160W/F606W flux-ratio map with intensity contours 
is displayed as an inset. 

Fig. 5.— Near-infrared images of IRAS 17245-3951: (a) F160W, (b) F222M, and (c) flux-ratio 
map of F160W/F606W, all on a logarithmic scale. The redder regions are darker on the flux-ratio 
map. The intensity range in (a) is 10"^^ to 3xl0~^ ergs cm~^ s~^ /xm arcsec"^ and in (b) is 
10~^^ to 3x10""^'^ ergs cm~^ s~^ /xm ~^ arcsec"^. 

Fig. 6.— Images of IRAS 16594-4656: (a) F160W, (b) F222M, (c) F606W, and (d) flux-ratio 
map of F160W/F606W, all on a logarithmic scale. The F606W image is included to show the faint 
nebulosity seen in the visible image, and it has the same orientation and scale as the F160W and 
F222N images. The intensity range in (a) is lO"-*^^ to 3xl0~^ ergs cm~^ /xm ~^ arcsec"^ and 
in (b) is 3xl0~^^ to 10~^ ergs cm~^ s~^ /itm ~^ arcsec~^. 

Fig. 7. — The 2 /im polarimctric results of the PPNs. The total intensity of the light through the 
three polarizers (Stokes I) is plotted by contours and the polarimetric vectors are plotted as lines 
whose lengths indicate the percentage polarization. The contour levels range from 98% to 0.1% 
of the peak values and are in equal intervals on a logarithmic scale. A scale for the percentage 
polarization is given for each object. The relative scales on the x and y axes are in arcsecs. 
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Fig. 8. — Polarimetric results for IRAS 17150—3224, displayed in different ranges of the percentage 
polarization: (a) 70% > p> 60%, (b) 60% > p> 50%, (c) 50% > p> 40%, (d) 40% > p> 30%, 
(e) 30% >p> 20%, (f) 20% >p> 10%, (g) p < 10%. These show a striking point symmetry. The 
intensity contours and x and y axis scales are as in Fig. 7. 



Fig. 9. — Polarimetric results for IRAS 17441—2411, displayed in different ranges of the percentage 
polarization: (a) p > 50%, (b) 50% > p > 40%, (c) 40% > p > 30%, (d) 30% > p > 20%, (e) 
20% > p > 10%, (f) p < 10%. A general point symmetry is seen. The intensity contours and x and 
y axis scales are as in Fig. 7. 



Fig. 10. — Polarimetric results for IRAS 17245—3951, displayed in different ranges of the percentage 
polarization: (a) p > 30%, (b) 30% > p> 20%, (c) 20% > p> 10%, (d) p < 10%. The intensity 
contours and x and y axis scales are as in Fig. 7. 



Fig. 11. — Polarimetric results for IRAS 16594—4656, displayed in different ranges of the percentage 
polarization: (a) 60% > p> 40%, (b) 40% >p> 30%, (c) 30% > p> 20%, (d) 20% > p> 10%, 
(e) p < 10%. The intensity contours and x and y axis scales are as in Fig. 7. 



Fig. 12.— 



Percentage polarization maps for IRAS 17150-3224 and IRAS 17441-2411. 



This figure "su.figl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0304400vl 



This figure "su.fig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0304400vl 



This figure "su.fig3a.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0304400vl 



This figure "su.fig3b.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0304400vl 



This figure "su.fig4.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0304400vl 



This figure "su.fig5.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0304400vl 



This figure "su.fig6.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0304400vl 



This figure "su.fig7.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0304400vl 



This figure "su.fig8.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0304400vl 



This figure "su.fig9.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0304400vl 



This figure "su.figl0.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0304400vl 



This figure "su.figll.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0304400vl 



This figure "su.figl2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0304400vl 



